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I.  Introduction 

This  report  presents  a  collision-induced  dissociation  (CID)  study  of  small  boron  oxide 
cations,  B„On,  motivated  by  a  need  for  more  accurate  and  reliable  structural  and  thermodynamic 
information  on  both  neutral  and  ionic  boron  oxides.  This  work  compliments  ongoing  ab  initio 
calculations  by  Page  and  co-workers,  and  provides  many  points  of  comparison  between 
experiment  and  theory.  In  addition,  the  results  of  our  CID  study  are  important  in  interpretation  of 
boron  oxide  cluster  chemistry  work  in  our  lab. 

In  combustion  of  boron  or  boron  rich  fuels  and  propellants,  boron  oxide  chemistry  plays 
two  important  roles.  Boron  particles  are  naturally  coated  with  a  passivating  layer  of  the  oxide, 
which  may  be  at  least  partly  converted  to  hydroxide  (i.e.  boric  acid)  depending  on  humidity  and 
temperature.  As  the  particles  heat  up  in  a  combustion  environment,  the  oxide  layer  retards  boron 
ignition  until  it  is  removed.  This  occurs  by  evaporation  at  high  temperatures,  but  the  potential 
exists  to  accelerate  ignition  by  chemical  volatilization  at  lower  temperatures.  Boron  oxide 
properties  are  also  important  in  achieving  energy  release  in  the  post-oxidation  chemistry.  In 
particular,  full  energy  release  only  occurs  if  the  gaseous  B„0„H|  oxidation  products  condense  to 
thermodynamically  stable  products  such  as  B203(?),  rather  than  remaining  as  (B„0„,H|)jj 
oligomers.  Cleariy,  the  chemistry  and  thermodynamics  of  small  B„0^i  species  are  important  in 
this  process. 

19970513  045 


WIC  QUALiry  INSPECTED  1 


k 


A  promising  strategy  to  enhance  boron  combustion  is  to  fluorinate  the  hydrocarbon 
components  of  a  boron-containing  propellant.  Decomposition  of  the  hydrocarbon  releases  HF, 
that  can  attack  the  boron  or  boron  oxide  particle  surfaces.  Since  boron  has  a  high  affinity  for 
fluorine,  and  F  is  isoelectronic  Avith  OH,  fluorine  will  tend  to  displace  OH  or  O  from  oxides, 
forming  stable  B„0„H|Ffc  species.  Since  fluorine  compounds  of  boron  are  more  volatile  than 
oxides  or  hydroxides,  HF  attack  should  tend  to  volatilize  the  oxide  layer.  Indeed  FBO  is  observed 
as  the  major  product  ofHF  attack  on  small  boron  oxide  cluster  ions  (1),  in  accord  with  modeling 
predictions  based  on  thermodynamic  considerations  (2,3). 

In  the  course  of  the  HF^oron  oxide  cluster  study,  we  found  substantial  uncertainties  in 
the  thermodynamics  reported  for  small  B„0„H|Fk  compounds.  Much  of  the  thermochemistry  is 
based  on  very  low-level  ab  initio  or  semi-empirical  quantum  chemistry  calculations,  that  are  not 
expected  to  be  very  reliable.  Recent  and  on-going  high  quality  ab  initio  calculations  by  Linder 
and  Page  (4,6)  and  Soto  (5)  are  improving  the  situation  dramatically,  however,  there  is  little 
experimental  data  available  to  test  the  accuracy  of  the  calculations.  The  results  reported  here 
largely  validate  the  accuracy  of  the  Linder  and  Page  calculations,  while  their  results  provide 
considerable  insight  into  our  experiments. 

One  major  puzzle  remains.  In  the  course  of  our  HF/B„0^  study  we  observed  a  reaction 
producing  FBOH.  From  the  energy  dependence  of  the  cross  section  for  this  reaction,  it  clearly  is 
exothermic  and  based  on  estimates  of  the  thermodynamics  for  B^O,^,  we  predicted  that  FBOH 
should  have  a  AHf  no  higher  than  -194  kcal/mole.  This  is  80  kcal/mole  more  stable  than  predicted 
by  the  calculations  of  Page  and  co-workers  (6)  and  Soto  (5),  and  this  discrepancy  is  far  outside 
the  usual  uncertainties  expected  for  calculations  of  that  quality.  Our  AHf  was  based  in  part  on 
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estimates  for  the  thermochemistry  of  the  larger  B„0^  species,  and  this  seemed  to  be  a  likely 
source  of  the  disagreement.  One  motivation  for  the  present  CID  study  was,  therefore,  to 
determine  the  stabilities  of  these  species  directly.  Despite  overall  good  agreement  between  our 
experiments  and  the  calculations,  this  particular  controversy  has  not  been  resolved. 

11.  Experimental  Method 

Apparatus.  The  cluster  beam  instrument  and  operating  procedures  used  for  these 
experiments  has  been  described  in  detail  elsewhere  (1,7,8).  Briefly,  boron  oxide  cluster  ions  are 
generated  by  12  keV  argon  atom  bombardment  (9)  of  a  film  of  vitreous  B2O3  maintained  near  its 
melting  point.  The  nascent  B^Oy^  cluster  ions  are  collected  by  a  radio-fi’equency  (rf)  octapole  ion 
guide  and  cooled  to  near  room  temperature  by  storage  in  a  labyrinthine  rf  trap  containing  -0.01 
Torr  of  helium  buffer  gas.  The  reactant  cluster  ion  size  and  composition  is  selected  using  a 
quadrupole  mass  filter  and  then  the  beam  is  injected  into  another  octapole  ion  guide  system  where 
scattering  is  carried  out.  The  octapole  sets  the  collision  energy  and  guides  the  ions  through  a 
collision  cell  filled  with  either  xenon  or  argon  to  a  typical  pressure  of  1  x  10'^  Torr.  For  these 
experiments  we  need  single  collision  conditions,  and  this  was  checked  by  measuring  cross  sections 
over  a  range  of  xenon/argon  pressures.  Fragment  ions  and  the  remaining  parent  ions  are  collected 
by  the  octapole,  mass  analyzed  by  a  second  quadrupole  mass  spectrometer,  and  counted. 

Sample  Preparation.  Isotopically  purified  (94. 1 1  At.%  of  ^®B)  boron  oxide  (  B2O3) 
powder  (Eagle-Pitcher)  was  sprinkled  on  a  stainless  steel  substrate,  then  heated  in  a  furnace  at 
650®C  for  about  four  hours  in  oxygen  environment  to  produce  a  vitreous  film.  Because  boron 
oxide  is  found  to  be  highly  hygroscopic,  the  sample  is  maintained  a  350-450‘’C  in  the  high 
vacuum  cluster  source  chamber. 
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III.  Results 


Reactant  Beam  Kinetic  Energy  Distribution,  Accurate  CE>  threshold  determinations 
require  a  narrow  and  well  defined  kinetic  energy  distribution  for  the  primary  beam.  This  is 
measured  by  retarding  potential  analysis  and  controlled  by  fine  adjustments  of  the  focussing  lens 
system.  A  typical  retarding  potential  curve  is  shown  in  Figure  1  correspond  to  a  beam  energy 
width  of  ~  0.240  eV.  The  beam  energy  distribution  is  obtained  by  fitting  this  retarding  curves  to  a 
3  parameter  asymmetric  Lorentzian  function; 

F(x)  = - ^ - 

{ l+[(x-Eo)/b]^  }exp(x-Eo)® 

Eo  represents  the  shift  in  average  beam  energy  relative  to  the  nominal  LAB  KE  of  the  ions,  the 
energy  width  is  determined  by  b,  and  c  is  allows  fitting  the  asymmetry  of  the  distribution  (usually 
quite  small). 

Reaction  Cross  Sections  and  Branching  Patterns.  Absolute  fi'agmentation  cross 
sections  were  measured  for  all  observable  (n<4,  m<5)  parent  ions  at  center-of-mass 
collision  energies  ranging  fi’om  0.8  to  10  eV.  Cross  sections  are  shown  for  all  significant 
fragmentation  channels  in  Figures  2a  through  2g.  Fragmentation  in  collision  with  both  argon  and 
xenon  was  studied  and  the  best  results  are  presented.  Xenon  generally  seems  to  be  a  more 
efficient  target  gas,  however,  argon  was  found  to  be  superior  for  several  of  the  smallest  cluster 
ions.  To  summarize  trends  with  changing  cluster  size.  Figure  3  reports  the  total  fi'agmentation 
cross  section  and  the  fi-agment  patterns  at  a  fixed  collision  energy  of  9.5  eV  —  well  above  the 
threshold  range  for  all  cluster  sizes. 

Extracting  CID  Thresholds.  To  obtain  quantitative  dissociation  thresholds  from  the 


4 


experimental  data  we  need  to  correct  for  the  collision  energy  spread  resulting  from  the  beam 
energy  distribution  and  the  thermal  motion  of  the  target  gas.  This  was  done  by  a  standard 
convolution  and  fitting  approach  that  has  been  described  in  detail  previously  (10).  Briefly,  the 
cross  sections  in  the  threshold  region  were  modeled  using  an  assumed  “true”  cross  section 
functional  form  that  has  been  widely  used  in  the  scattering  community: 


A(E^vaa  -Eo)“ 

®C^avail)  ~ 


Eavail 

where  ^4  is  a  normalization  factor,  Eq  is  the  dissociation  threshold,  /i  is  an  adjustable  parameter 
that  varies  the  curvature  of  the  function,  physically  related  to  the  energy  transfer  efficiency,  and 
^avcdi  is  the  total  energy  available  to  drive  fragmentation. 

This  trial  function  is  convoluted  (through  a  Monte  Carlo  simulation)  with  the  kinetic 
energy  distribution  of  the  primary  cluster  ions,  the  translational  energy  distribution  of  the  target 
gas,  and  the  distribution  of  vibrational  and  rotational  energy  of  cluster  ions.  The  vibrational 
energy  has  been  calculated  assuming  that  the  clusters  are  at  400  K,  using  vibrational  frequencies 
from  ab  initio  calculations  (6).  The  parameters  n  and  Eq  are  optimized  until  a  best  fit  is  obtained. 
Fitting  was  attempted  only  for  the  lowest  energy  fragmentation  channel  for  each  reactant  cluster 
ion.  In  principle,  additional  thermochemical  information  can  be  extracted  from  the  higher  energy 
dissociation  channels,  however,  this  analysis  is  complicated  by  a  difficult-to-estimate  kinetic  shift 
factor. 


The  best  fits  for  the  lowest  energy  fragmentation  channels  are  shown  in  Figures  2a  -  2g. 
The  curves  labeled  “Best  fit”  are  best  free  fits  to  the  data,  and  the  extracted  Eq  values  are  the 
experimental  estimates  for  dissociation  energies.  We  also  plot  "Fit  to  Theory"  curves.  These  are 
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fits  based  on  Eq  values  calculated  by  Linder  and  Page  (6)  with  only  the  n  parameter  adjusted  in  an 
attempt  to  fit  the  data.  These  are  plotted  to  show  the  cases  for  which  the  experiment  and  theory 
are  or  are  not  in  agreement.  The  fit  thresholds  (Eq)  and  n  parameters  for  the  lowest-energy 
dissociation  channel  of  each  cluster  for  both  types  of  fittings  are  given  in  Table  I,  along  with  the 
calculated  dissociation  thresholds.  The  cluster  stability,  represented  by  the  lowest-energy 
dissociation  thresholds,  is  plotted  as  a  function  of  size  in  Figure  3  (scatter  plot,  right-hand  scale). 
In  addition,  stability  is  plotted  as  a  function  of  boron-to-oxygen  ratio  in  Figure  4. 

IV.  Discussion 

Fragmentation  Branching  Ratios.  In  fi-agmentation  reactions,  two  factors  can  contol  the 
branching  between  different  possible  fragment  channels.  In  some  cases,  one  might  expect  that  the 
geometrical  struture  of  the  parent  cluster  ion  would  be  reflected  in  the  fragments,  i.e.  that 
fi’agmentation  might  occur  by  simple  bond  rupture,  without  rearrangement.  This  is  particularly 
likely  in  high  energy  CID,  where  the  fragmentation  time  scale  is  short.  In  low  energy  CID,  such 
as  our  near-threshold  work,  the  fi'agmentation  time  scale  is  long  enough  that  rearrangement  can 
occur  prior  to,  or  during,  decomposition.  In  this  case,  the  fi'agment  branching  is  influenced 
strongly  by  the  thermochemistry  of  the  possible  products,  with  the  most  stable  products 
dominating. 

Based  on  the  available  thermochemistry  in  the  literature  (1,4,5, 1 1)  and  the  on-going  ab 
initio  calculations  by  Linder  and  Page,  it  appears  that  product  thermochemistry  is  the  dominant 
factor  for  all  but  one  of  the  boron  oxide  cluster  ions  we  have  examined.  For  BO^,  BjO^,  B2O2*, 
and  B303^  —  all  clusters  where  the  number  of  B  atoms  is  equal  to  or  greater  than  the  number  of  O 
atoms  in  the  parent  cluster,  the  dominant  fragmentation  pathway  is  loss  of  B^.  This  appears  to 
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reflect  two  factors:  the  low  ionization  potential  of  the  boron  atom,  and  the  relatively  high  stability 
of  fi'agments  with  stoichiometry. 

For  the  8203^  parent  cluster,  the  only  important  fragmentation  channel  is  loss  of  O, 
yielding  6202^-  Dissociation  to  BO^  +  BO2,  only  0.5  eV  higher  in  energy,  accounts  for  only  a 
few  percent  of  the  products,  even  at  energies  well  above  threshold.  8364^  which  also  prefers  to 
decompose  to  8202^  (+  BO2),  produces  a  substantial  branching  to  the  nearly  isoenergetic 
fragment  pairs  8263^  +  80  and  80'’  +  B2O3  that  lie  -0.3  higher  in  energy. 

The  one  exception  to  control  by  fragment  thermochemistry  is  802^.  The  dominant 
decomposition  charmel  is  to  80^  +  O,  even  though  the  8^  +  O2  channel  is  - 1 .6  eV  lower  in 
energy.  This  presumably  reflects  either  a  barrier  or  dynamical  bottleneck  that  inhibits  passage 
from  the  parent  080  structure  to  a  transition  where  O2  elimination  can  occur.  Some  8^  is 
observed  with  approximately  the  same  appearance  energy  as  the  main  80^  channel.  The 
competition  between  these  channels  makes  analysis  of  the  80^  +  O  threshold  somewhat 
ambiguous. 

Given  that  thermochemistry  appears  to  control  the  fragment  branching,  some  energetic 
insight  can  be  inferred  directly  from  the  branching  patterns.  In  several  cases  we  observe  pairs  of 
channels  that  differ  only  in  which  fragment  carries  the  charge: 

NT  -*  A"  +  8 

-*  A  +  B". 

As  the  fragments  separate,  we  expect  that  the  charge  will  largely  end  up  on  the  fragment  with  the 
lower  ionization  potential  (IP).  If  the  IPs  are  similar  both  pairs  of  channels  will  be  observed,  but 
if  the  difference  is  large,  we  expect  significant  signal  only  for  the  lower  energy  fragment  pair. 
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Based  on  this  assumption,  we  can  infer  relationships  between  IPs  for  a  number  of  B„Om  species 
and  these  are  summarized  in  Table  I.  In  many  cases  these  are  not  surprising  and  agree  with 
literature  thermochemistry.  In  others,  the  literature  values  are  inconsistent  with  our  results  or  no 
literature  values  are  available.  The  on-going  calculations  of  Linder  and  Page  appear  to  be 
consistent  with  our  findings  in  all  cases. 

One  notable  example  is  B2O3.  The  existing  literature  gives  an  IP  for  B2O3  of  13.56  eV 
(1 1).  Our  fi’agmentation  branching  for  B304^  indicates  that  IP(B203)  <  IP(BO),  given  as  13  .0  eV 
in  the  literature.  Linder  and  Page  have  calculated  IP(BO)  =  13.1  eV  and  ^(8203)  =  12.9  eV, 
consistent  with  our  findings. 

Fragmentation  Threshold  Energies  and  Relation  to  Theory. 

These  CID  experiments  show  that  B02^  is  the  most  stable  and  830^  the  least  stable 
cluster  studied.  The  stability  of  the  boron  oxide  clusters  plotted  as  a  function  of  the  cluster  size 
(Figure  3)  exhibits  an  oscillating  pattern  typical  for  most  of  size  dependence  cluster  studies 
(8,10).  Furthermore,  the  stability  clearly  is  strongly  anticorrelated  to  the  B:0  ratio  in  the  cluster 
(Figure  4).  As  expected  maximum  stability  is  promoted  by  increasing  the  possibility  for  BO 
bonding. 

Table  II  gives  the  best  fit  dissociation  threshold  extracted  as  discussed  above,  along  with 
the  associated  value  of  the  "n"  parameter.  Also  given  are  the  dissociation  energies  calculated  by 
Linder  and  Page.  Note  that  the  experimental  uncertainties  are  abnormally  large  for  these  oxide 
cluster  ions.  This  results  fi'om  the  fact  that  the  fragmentation  cross  sections  are  unusually  small 
(at  maximum,  only  a  few  percent  of  the  collision  cross  sections)  and  rise  rather  slowly  fi'om 
threshold.  This  suggests  that  collisional  energy  transfer  is  particularly  inefficient  for  these 
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clusters.  The  result  is  a  larger-than-usual  degree  of  arbitrariness  in  picking  the  best  fit  in  the  near¬ 
threshold  energy  range  because  the  signal  rise  very  slowly  out  of  the  background.  To  give  a 
clearer  idea  of  which  experimental  dissociation  threshold  (Eq)  values  agree/disagree  with  the 
calculated  dissociation  energies,  figures  2a-2g  plot  the  best  fits  obtainable  using  the  calculated  Eq 
values. 

For  BO^  there  is  an  apparent  discrepancy  between  the  experimental  and  calculated  Eo, 
however,  as  figure  2a  shows,  the  "best"  and  "theory"  fits  are  rather  similar.  BO^  has  an  especially 
small  CID  cross  section,  and  the  data  are  not  good  enough  to  justify  concluding  that  the 
discrepancy  is  real. 

For  B20^  and  some  of  the  larger  clusters  there  is  an  additional  complication.  Linder  and 
Page  calculate  that  the  BOB^  isomer  is  most  stable,  but  that  cyclic  and  BBO^  isomers  lie  only  0.2 
eV  and  0.3  eV  higher  in  energy.  For  this  situation,  it  is  not  unlikely  that  our  cluster  beam  contains 
some  admixture  of  these  higher  energy  isomers,  which  would  give  a  threshold  of  ~  1 .6  eV,  in  good 
agreement  with  the  experiment.  Because  CID  is  inherently  sensitive  to  the  least  stable  species 
present,  only  a  few  percent  of  a  high  energy  isomer  is  needed  to  skew  the  threshold.  This  is  NOT 
the  case  for  reactive  scattering,  where  usually  we  are  looking  a  far  more  efficient  process  than 
CID. 

This  example  raises  the  issue  of  the  isomer  distributions  found  in  our  cluster  ion  beams. 
Our  cluster  ions  are  produced  by  laser  ablation,  or  in  this  case,  particle  sputtering.  Both  methods 
produce  a  distribution  of  cluster  sizes  with  very  high  internal  temperatures.  These  hot  cluster  ions 
initially  cool  by  evaporation,  then  are  actively  cooled  by  coUisional  quenching  in  helium  buffer  gas 
over  a  ~  1  second  period.  Based  on  previous  tests,  we  are  pretty  confident  that  the  clusters  have 
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internal  temperatures  near  room  temperature.  As  the  clusters  cool  down,  we  normally  would 
expect  that  nearly  all  the  clusters  will  anneal  to  the  lowest  energy  isomer,  since  the  density  of 
states  is  highest  in  the  deepest  well.  Indeed,  in  our  previous  studies  of  boron  and  carbon  cluster 
ions,  the  evidence  (agreement  with  theory,  fragment  distributions,  and  reactivity)  indicates  that 
with  the  exception  of  a  few  particular  carbon  clusters,  only  the  most  stable  isomer  is  found.  There 
are  two  cases  where  a  substantial  concentration  of  a  higher  energy  isomer  may  be  present.  In  the 
case  of  C7.9,  the  linear  and  cyclic  isomers  are  nearly  isoenergetic,  so  that  densities  of  states  for 
each  form  are  comparable. 

The  other  likely  case  is  if  the  activation  barriers  separating  the  potential  minima 
corresponding  to  the  two  isomers  are  large  compared  to  the  energy  difference  between  them.  In 
this  case  the  barrier  crossing  rate  will  decrease  rapidly  early  in  the  cooling  process,  thus  freezing 
in  an  isomer  distribution  corresponding  to  a  temperature  far  above  the  final  cluster  temperature. 
For  some  of  these  boron  oxide  clusters,  this  may  indeed  be  the  case.  The  clusters  are  generally 
singly  connected  (unlike  boron  for  example)  so  that  the  bond  order  reduction  in  fragmentation  is 
similar  to  that  in  the  type  of  transition  state  that  seems  likely  for  the  isomerization.  This  issue  is 
being  examined  by  Linder  and  Page,  who  are  calculating  isomerization  activation  barriers  for 
some  of  the  clusters  where  multiple  isomers  seem  likely.  If  this  does  turn  out  to  be  a  problem  in 
the  experiments,  we  can  attempt  to  remove  the  higher  energy  isomers  by  titrating  them  out  during 
the  cooling  process  using  a  suitable  reactant  -possibly  water.  This  approach  has  been  used 
successfully  for  carbon  cluster  ions. 

For  B02^  (unlikely  to  have  any  stable  isomers  other  than  OBO^)  our  best  fit  Eq  is 
somewhat  lower  than  the  calculated  Eg,  but  as  in  the  BO^  case,  this  cluster  has  a  small 
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fragmentation  cross  section  and  the  experimental  uncertainty  is  therefore  high. 

8202^  is  really  the  first  cluster  where  there  clearly  is  a  difference  between  the  data  and 
calculation.  The  experimental  signal/noise  is  better  here,  and  the  "fit  to  theory"  does  not 
adequately  reproduce  the  threshold  behavior.  It  is  possible  that  this  reflects  a  problem  in  the 
calculations,  however,  8262^  is  another  case  where  a  low  lying  isomer  might  skew  the 
experiments.  Linder  and  Page  predict  that  there  is  a  cyclic  isomer  lying  -0.8  eV  above  the  most 
stable  0808*  geometry.  This  would  give  a  threshold  of  3  .3  eV,  consistent  with  the  measured 
Eo- 

82O3*  is  a  case  where  isomers  other  than  the  ground  state  08080*  are  not  expected  to 
be  low  enough  lying  to  be  significantly  populated.  For  this  cluster  the  experimental  and  calculated 
Eq  values  are  in  nearly  perfect  agreement. 

83O3*  is  another  case  where  it  is  not  possible  to  fit  the  experimental  data  with  the  Eq  value 
calculated  by  Linder  and  Page.  The  experimental  Eq  is  over  1  eV  lower  than  the  calculation.  As 
with  all  the  oxides  with  8:0  ratio  equal  to  or  greater  than  unity,  several  low  lying  isomers  are 
likely.  The  calculations  find  the  088080*  structure  to  be  most  stable,  but  080808*  and  a 
structure  with  a  4-membered  ring  are  both  calculated  to  be  within  - 1  eV  of  the  ground  state.  A 
smaU  admixture  of  one  or  both  of  these  might  explain  the  near-threshold  discrepancy  shown  in 
figure  2f 

For  83O4*  both  intuition  and  the  Linder-Page  calculations  suggest  a  0808080* 
structure.  Calculations  are  on-going,  but  to  date  no  other  low-lying  isomers  have  been  found. 
Nonetheless  there  is  a  substantial  discrepancy  between  the  experimental  Eq  (4.25  eV)  and  that 
derived  from  the  calculations  (6.4  eV).  At  this  point,  the  origin  of  the  disagreement  is  not  clear. 
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Implications  for  FBOH  thermochemistry.  Another  unresolved  discrepancy  is  the  heat  of 
formation  of  FBOH.  In  our  studies  of  boron  oxide  cluster  ion  reactions  with  HF,  one  of  the 
major  product  channels  observed  was  elimination  of  FBOH: 

B„0„^  +  HF  -^B„.iO„/  +  FBOH 
-♦B„.iO„.i+FBOH" 

Production  of  neutral  FBOH  is  inferred  from  the  B„.iO„.i^  product,  for  which  FBOH  is  the 
lowest  energy  neutral  partner.  In  reaction  with  B304^,  both  these  channels  are  observed  with 
nearly  equal  intensities  and  with  substantial  cross  sections  at  low  collision  energies,  indicating  that 
both  are  exothermic.  The  observation  that  the  reaction  producing  B„.iO„.i^  (+  FBOH)  is 
exoergic  allows  us  to  put  an  upper  limit  on  AH{(FBOH).  Taking  Linder  and  Page's  calculated 
AHf  values  for  B304^  (-  44.2  kcal/mole)  and  B203^  (100.8  kcal/mole),  this  requires  that 
AHjCFBOH)  <210  kcal/mole.  This  compares  poorly  with  the  ab  initio  AH£(FBOH)  =  —113 
kcal/mole  calculated  by  Soto  (5)  and  Page  (4).  Our  CID  results  find  B304'^  to  be  ~49  kcal/mole 
less  stable  than  calculated  by  Linder  and  Page.  If  we  use  this  to  estimate  AHf(B304^)  =  ~  5 
kcal/mole,  this  still  requires  that  AHf(FBOH)  <  ~  160  kcal/mole.  In  either  case,  the  discrepancy 
with  the  ab  initio  value  is  way  outside  the  range  of  errors  normally  associated  with  those 
calculations. 

The  idea  that  FBOH  is  considerably  more  stable  than  the  ab  initio  calculations  suggest  is 
supported  by  the  observation  that  the  (Bn.iO^.!"^  +  FBOH)  and  (B„.iO„.i  +  FBOIT)  product 
channels  have  nearly  equal  intensity  for  reaction  of  B304^  with  HF.  This  suggests  that  IP(FBOH) 
=  IP(B203)  =  12.9  eV  (6).  The  difference  between  ab  initio  AHf  values  for  FBOH  and  for 
FBOH^  (the  latter  being  consistent  with  our  results)  gives  IP(FBOH)  =  ~7.5  eV.  If  this  were 
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correct,  it  seems  unlikely  that  we  would  see  significant  signal  for  the  +  FBOH  channel 

because  the  charge  should  overwhelmingly  migrate  to  FBOH  as  the  products  separate. 
Conclusions 

Our  CID  results  help  clarify  some  of  the  confusion  in  the  thermochemical  literature  for 
small  boron  oxide  oligamers.  The  results  are  largely  consistent  with  the  new  ah  initio  calculations 
of  Linder  and  Page,  though  there  are  some  unresolved  issues  regarding  isomer  distributions  that 
make  direct  comparisons  ambiguous  for  some  cluster  sizes.  Significant  (and  related) 
discrepancies  between  theory  and  experiment  exist  for  8304^  and  for  FBOH.  As  the  calculations 
are  completed  for  6304%  it  may  be  possible  to  resolve  these  disagreements. 
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Table  L 


Comparisons  of  the  Ionization  Potentials  of  Boron  oxide  Cluster  Ions. 


parent  oxide 
cluster 

IP  relations  inferred 

literature  IP 
relations^ 

Linder  and  Page  IP 
relations 

IP(BO) » IP(B) 

IP(BO)  =  13.0eV 
IP(B)  =  8.298eV 
consists 

IP(BO)=13.1eV 
IP(B)  =  8.30eV 

IPCBOj)  >  IP(B) 

IP(B02)  =  13.5eV 
consistent 

IP(B02)  =  13.66eV 

B,Of 

IP(BO)«  IP(B02) 

IP(BO)  «  IP(B02) 

IP(B  A)  » IP(B) 

IP(B203)  =  13.56eV 
consistent 

IP(B203)  =  13.15eV 

IP(BO)  >  IPCBjOj) 

IP(B202)  =  13.58eV 
inconsistent 

IP(B202)  =  10.37eV 

B3O/ 

IP(BO)  >  IPOBjOj) 
IP(B02)  >  IP(B202) 

inconsistent 

inconsistent 

^  Lias,  S.  G.,  Bartmess,  J.  E.,  Liebman,  J.  F.,  Holmes,  J.  L.,  Levin,  R.  D.,  Mallard, 
W.  G.  J.  Phys.  Chem.  Ref.  Data,  1988,  vol.  17,  Suppl.l. 
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Table  n.  Best  Fit  parameters  (n)  and  Best  Fit  Dissociation  Energy  Thresholds 
(Eq)  of  Boron  Oxide  Cluster  Ions. 


dissociation 

best  fit 

^t  to  theory 

channel 

n 

Eo 

n 

Eo 

BO^  -  B^  +  O 

1.90 

2.75  ±0.6 

1.35 

3.60 

BjO^  -  B^  +  BO 

1.80 

1.30  ±0.5 

1.40 

1.90 

BOj^  -  BO^  +  O 

1.60 

4.75  ±0.6 

1.05 

5.60 

BjOj^  -  B^  +  BO2 

1.30 

3.15  ±0.2 

1.05 

4.10 

BjOj^  -  +  0 

1.64 

5.10  ±0.2 

1.70 

5.09 

BjOj*  -  B*  +  BjOj 

2.05 

2.00  ±0.2 

1.05 

3.2 

2.10 

4.25  ±0.5 

1.1 

6.4 

*  Eq  values  from  reference  6. 
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Figure  Captions 

1.  A  typical  retarding  potential  curve. 

2.  Cross  sections  for  all  significant  fragmentation  channels  along  with  the  best  fits  for  the  lowest 
energy  dissociation  pathways  of  boron  oxide  cluster  cations. 

3.  The  total  fragmentaion  cross  section  and  the  fragmentation  patterns  at  a  fixed  collision  energy 
of  10  eV.  The  cluster  stability  as  a  function  of  size  is  shown  by  the  scatter  plot  (right-hand  scale). 

4.  The  cluster  cations’s  stability  as  a  function  of  boron-to-oxygen  ratio. 
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